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Abstract 



The process e"'"e~ — > 77(7) is studied using data collected by the OPAL detector at LEP 
between the years 1997 and 2000. The data set corresponds to an integrated luminosity of 
672.3 pb~^ at centre-of-mass energies lying between 181 GeV and 209 GeV. Total and differential 
cross-sections are determined and found to be in good agreement with the predictions of QED. 
Fits to the observed angular distributions are used to set limits on parameters from several 
models of physics beyond the Standard Model such as cut-off parameters, contact interactions 
of the type e^e~77, gravity in extra spatial dimensions and excited electrons. In events with 
three photons in the final state the mass spectrum of photon pairs is investigated. No narrow 
resonance X — >■ 77 is found and limits are placed on the product of the X7 production cross- 
section and branching ratio. 
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1 Introduction 



This paper presents a study of the process e+e~ 77(7), where the brackets indicate a 
possible third particle that might escape along the beam direction. Data with a total integrated 
luminosity of 672.3 pb~^ collected in the years 1997 to 2000 with the OPAL detector at LEP 
at centre-of-mass energies between 181 GeV and 209 GeV are used. This process is one of the 
few reactions which are dominated by QED even at these energies. The high statistics sample 
considered in this paper allows precision tests of QED and searches for new particles such as 
excited electrons and photonically decaying resonances. 

This process has been studied previously at LEP at lower energies [|l],|^,0 and by other experi- 
ments 1^, |, 0, H, 1^, [ly, |n|] . The increased data sample since our last pubhcation leads to a 
much better understanding of the detector and the systematic errors. In addition the selection 
efficiency has been increased. Therefore this paper includes a re-analysis of previously published 
data at 183 GeV and 189 GeV and supersedes those results. Total and differential cross- 
sections are determined and compared with the predictions of QED. 



There are several models which predict deviations from the QED cross-section introducing cut- 
off parameters |12|, e^e~77 contact interactions |1^, the exchange of excited electrons |14, 15 



or the exchange of Kaluza-Klein gravitons in models with extra dimensions [O] . The measured 



differential cross-sections are used to place constraints on the parameters of these models. 

The search for a resonance X which is produced along with a photon and decays via X ^ 77 is 
especially interesting in the light of models with fermiophobic Higgs bosons [|l^] or hypercharge 
axions [18|. Because of the higher kinematic reach the three-photon final state is complementary 
to ZX production. Events with three observed photons are used to search for such a resonance. 
Similar searches were performed previously at the Z peak |jl9| and at lower LEP2 energies 



In the next section a short description of the relevant detector components is given. This is 
followed by an overview of the data and Monte Carlo samples used in the analysis. Section ^ 
contains an introduction to the theoretical models. The event selection is described in Section 
the systematic errors are summarised in Section |^ and the results are given in Section |^. 



2 The OPAL detector 



A detailed description of the OPAL detector can be found elsewhere so only those 

components most relevant to this analysis are discussed here. 

The most important subdetector for this analysis is the electromagnetic calorimeter (ECAL). 
The barrel regionQ (| cos^| < 0.82) consists of 9440 lead-glass blocks in a quasi-pointing geo- 
metry each with a cross-section of about 10 x 10 cm at an inner radius of 2.45 m. The two 

""^OPAL uses a right-handed coordinate system in which the z axis is along the electron beam direction, x 
points to the centre of LEP such that the y axis points approximately upwards. The polar angle, 6, is measured 
with respect to the z axis and the azimuthal angle, (j), with respect to the x axis. 
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endcaps (0.81 < | cos6\ < 0.98) consist of 1132 blocks each, aligned parallel to the beam axis. 
The inner faces are placed at a z position of about ±2.3 m. For beam-energy photons, the 
spatial resolution is about 11 mm, corresponding to an uncertainty of 0.2° in 6, and the energy 
resolution is about 2% in the barrel and 3% - 5% in the endcaps, depending on the polar angle. 

The ECAL surrounds the tracking chambers. The large volume jet chamber (CJ) has an outer 
diameter of 3.7 m and is about 4 m long. It is segmented into 24 sectors each containing 159 
axial sense wires. The vertex chamber (CV) has an outer diameter of 47 cm and is 1 m long. 
It is segmented into 36 sectors each containing 12 axial wires in the inner region and an outer 
region of 6 stereo wires. CV and CJ are inside a common pressure vessel and separated by a 
foil and a carbon fibre tube. A silicon micro-vertex detector (SI) is located between CV and 
the beam pipe, at a radius of about 7 cm. This subdetector consists of two layers of double 
sided [z and sensitive) ladders. The inner (outer) layer covers | cos6\ < 0.93 (0.89) and has 
a (j) acceptance of 97.8% (99.6%). 

Outside the ECAL the instrumented return yoke serves as a hadronic calorimeter. The out- 
ermost subdetector is the muon system which consists of up to 4 layers of drift chambers. 
The luminosity is measured using small-angle Bhabha events collected in the silicon-tungsten 
luminometer | |23| ] . 



3 Data sample and Monte Carlo simulation 

Table |l| lists the data samples analysed in this paper. They are from the last four years of LEP 
running with a total of 5235 selected events taken at centre-of-mass energies between 181 GeV 
and 209 GeV. As an example, the last e"'"e~ — > 77(7) event, which was recorded three minutes 
before the final shut down of LEP, is shown in Figure |l]. Two photons, back-to-back in 0, are 
observed in the detector. One photon has converted between CV and CJ. From the acoUinearity 
a third photon is deduced to have escaped along the beam direction. 

Various Monte Carlo samples are used to study the selection efficiency and expected back- 



year 


V^nom 


V^rangc [GeV] 


[GeV] 


luminosity [pb ^] 


events 


1997 


183 




180.8- 


-184.0 


182.68 


55.57±0.15±0.19 


538 


1998 


189 




188.6- 


-189.0 


188.63 


181.07±0.16±0.36 


1531 


1999 


192 




191.4- 


-192.0 


191.59 


29.03±0.06±0.07 


258 


1999 


196 




195.2- 


-196.0 


195.53 


75.92±0.10±0.16 


616 


99/00 


200 




199.4- 


-200.2 


199.52 


78.20±0.11±0.17 


554 


99/00 


202 




201.4- 


-202.5 


201.63 


36.78±0.07±0.08 


281 


2000 


205 




202.5- 


■205.5 


204.88 


79.22±0.11±0.17 


566 


2000 


207 




205.5- 


■209.0 


206.56 


136.49±0.14±0.30 


891 



Table 1: Data used in this paper. The year of data taking, the nominal centre-of-mass energy, 
the energy range, the luminosity-weighted mean centre-of-mass energy, the integrated luminos- 
ity with its statistical and systematic error and the number of selected events are shown. 
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ground contributions. The signal e+e" 77(7) events are generated using the RADCOR ["M 
generator which rehes on a full 0{a^) calculation taking the electron mass into account. Events 
with four observed photons can be simulated using FGAM which contains a lowest-order 
relativistic calculation of four-photon production. Bhabha events are simulated using BH- 
WIDE [g§ (e+e-(7)) and TEEGG 

NUNUGPV [H are used. Tau pairs are simulated using KORALZ [|8 
is used for hadronic events. All samples are processed through the OPAL detector simulation 
program 



events KORALZ @ and 
and PYTHIA |30|,^ 



32 1 and reconstructed in the same way as the data. 



4 Cross-section for the process e+e — > 77(7) 
4.1 QED Born cross-section 

Up to the highest LEP energies the process e"'"e~ — >■ 77(7) can be described by QED. The 
lowest order cross-section in the relativistic limit is given by: 

( da \ 2tx o? 1 + cos^^^ ^ 



d cos 6 I ^ s 1 — cos^ 9 

^ / Born 



where ^/s is the centre-of-mass energy and a is the fine-structure constant at zero momentum 



transfer. The non-relativistic cross-section formula is given in [Q. Since the final-state particles 
are identical, the polar angle 6 is defined such that cos^ > 0. At Born level the definition of 9 
is unambiguous since the two photons are back-to-back. The definition at higher orders used 
in this paper is given in Section [4.2|. 



Weak interactions contribute only via loop diagrams. At the W-pair threshold weak corrections 
of up to 1.2% are expected for cos6' = At the energies considered in this analysis the 

corrections are smaller, e.g. for a centre-of-mass energy of 200 GeV they are less than 0.2% for 
all angles, and will be neglected. 



4.2 Radiative corrections 



In the presence of higher order effects the measured angular distribution depends not only on one 
angle 6 but on the angles of all produced photons. To enable a comparison of the measurement 
with Equation |I| an event angle must be defined to substitute 6. Various definitions of the 
event angle are possible, each leading to a different measured angular distribution. The ratio of 
this physical distribution to the distribution from the Born-level prediction is called a radiative 
correction TZ. The event angle 6* used in this analysis is defined by: 




where 61 and 62 are the angles of the two highest-energy photons. At Born level cos^ = cos 6'*. 
The angle 6* is equivalent to the angle in the centre-of-mass system of the two highest-energy 
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photons unless a third photon is produced away from the beam direction. With this definition, 
deviations of TZ from unity are relatively small and uniform as determined from an 0{a^) Monte 
Carlo 1^ without detector simulation. Using high statistics samples of 10^ events, generated 
at centre-of-mass energies of 189 GeV, 200 GeV and 206 GeV, the correction TZ is found to 
depend only weakly on the centre-of-mass energy. Therefore the y/s dependence is neglected. 
The average correction for cos^^* < 0.93 is TZ = 1.0448; the angular dependence can be found 
in Table |^. For QED processes the effects due to the next order can be assumed to be about 
10% of the corrections at this order. Since no fourth order Monte Carlo generator is available 
and weak contributions are neglected, a systematic error of 0.01 is assumed for the radiative 
correction. The ratio TZ is used to correct the measured cross-sections presented in this paper 
to the Born level. 



4.3 Alternative models 

Various models predict deviations from the QED expectation^. The simplest ansatz is a 
short-range exponential deviation from the Coulomb field parameterised by cut-off parame- 
ters A± [jl2|. This leads to a differential cross-section of the form 



do- \ ( da \ , a^Tis ,^ o ^, 

±-— (l + cos^^^). (3) 



y d cos d ) \d cos 6 J g^^^ 

New effects can also be introduced in effective Lagrangian theory [|l^. Here dimension-6 terms 
lead to anomalous ee7 couplings. The resulting deviations in the differential cross-section are 
similar in form to those given in Equation |^, but with a slightly different definition of the 
parameter: Ag = |-A^. Dimension 7 and 8 Lagrangians introduce ee77 contact interactions 
and result in an angle-independent term added to the Born cross-section: 

( da \ ( do \ 1 

(4) 



\dcos9 J j^, \dcos9j^^^^ 16 A' 

The associated parameters are given by A7 = A' and Ag = mgA'^ for dimension 7 and dimen- 
sion 8 couplings, respectively. The subscript refers to the dimension of the Lagrangian. 

Instead of an ordinary electron, an excited electron e* with mass Mg* could be exchanged in 
the t-channel |I^, p!5|] . In the most general case e*e7 couplings would lead to a large anomalous 



magnetic moment of the electron |^,|3^. This effect can be avoided by a chiral magnetic 
coupling of the form 



r . = 6*0^"" 

2 A 



cl + h.c. , (5) 



where r are the Pauli matrices and Y is the hypercharge. The parameters of the model are 
the compositeness scale A and the weight factors / and /' associated to the gauge fields W 
and B with Standard Model couplings g and g' . For the process e+e" — > 77(7), the following 
cross-section results 



All cross-sections given here are to lowest order. 
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da \ I da \ 

^ d cos 6* y \d cos 9 J ^^^^ 

2 -f 4 

2 A4 ^ 



(6) 



with = — 1(/ + /'); = ~f (1 ~ cos6') and = — 1(1 + cos^^). Effects vanish in the case 
of / = -/'. 

Theories of quantum gravity in extra spatial dimensions could solve the hierarchy problem 
because gravitons would be allowed to travel in more than 3+1 space-time dimensions |]39[ . 
While in these models the Planck mass Md in D = n + 4 dimensions is chosen to be of 
electroweak scale the usual Planck mass Mpi in four dimensions would be 

= i?"M]^+2 , (7) 

where R is the compactification radius of the additional dimensions. Since gravitons couple 
to the energy-momentum tensor, their interaction with photons is as weak as with fermions. 
However, the huge number of Kaluza-Klein excitation modes in the extra dimensions may give 
rise to observable effects. These effects depend on the scale Ms(~ Md) which may be as low 
as C(TeV). Model dependencies are absorbed in the parameter A which is expected to be of 
For this analysis it is assumed that A = ±1. The expected differential cross-section is 
given by [[T6| : 

\dcos9j \dcos9j^^^^ Stt M| 



5 Event selection 



Multi-photon events have a very clear experimental signature. They have large energy deposits 
in the electromagnetic calorimeter and small missing transverse momentum. They share these 
properties with Bhabha events (e"'"e~ e"'"e~(7)), which represent the most important back- 
ground. The first stage of the event selection is a preselection which requires events to have 
photon candidates and a low multiplicity of tracks and clusters. Cuts are then applied to reject 
cosmic-ray background. This is followed by a kinematic selection based on the signature of 
the energy deposit in the ECAL. The last step is a neutral event selection designed to reject 
Bhabha events using information from the tracking chambers. After application of all selection 
criteria the total background is reduced to an almost negligible level. 



5.1 Preselection 



A photon candidate is defined here as an ECAL cluster with an energy of at least 1 GeV and 
a polar angle satisfying | cos^| < 0.97. The cluster must consist of at least two ECAL blocks 
to ensure a good determination of the photon angle. 
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The preselection is very loose, requiring at least two photon candidates and a total ECAL 
energy in the event of at least O.ly/s. High multiplicity events are rejected by restricting the 
sum of the number of tracks and ECAL clusters to < 17. 

The two highest-energy photons must be within | cos^| < 0.93. This avoids a region where the 
detector material is not sufficiently well modelled to ensure a good description of the conversion 
probability and the angular reconstruction. The search for a resonance X is based on three- 
photon events and the invariant mass of X is calculated from the photon angles. Only events 
in which all three photons are in the region | cos^| < 0.93 are used for this search. 

5.2 Cosmic- ray background rejection 

A cosmic-ray particle can create signals in the outer detectors without producing a reconstructed 
track in the central tracking chambers. These particles do not necessarily pass close to the beam 
axis. Since the hadronic and electromagnetic calorimeters have different radii, the resulting hits 
in the two detectors occur separated in azimuth. Background of this type is suppressed using 
information from the muon chambers and the hadronic calorimeter. Events are rejected if there 
are three or more track segments reconstructed in the muon chambers. Events are also rejected 
if the highest-energy hadronic cluster is separated from each photon candidate by at least 10° 
in (p and has an energy of more than 20% of the summed energy of the photon candidates. In 
the case of one or two muon track segments the latter cut is tightened to 10%. Additionally, 
events are rejected if there is a track with a momentum of more than 10 GeV separated from 
each of the photon candidates by more than 10° in 0. 

5.3 Kinematic event selection 

Events are selected if they have large total energy and small missing transverse momentum. 
To improve the resolution on these quantities the information on the angle of the cluster is 

used instead of the energy information where possible. To facilitate this, the events are divided 
into four classes (/, //, III and IV) according to the number of photon candidates N^, the 
acoUinearity ^acoi = 180° — and the aplanarity Kaplan = 360° — {au + ctia + 0:23), where aij 



Cacol < 10° 


Uoi > 10° 


/ 


II 




III 


Kaplan < 0.1° 


= 3 


IV 


Caplan > 0.1° 


A^ > 4 



Table 2: Definition of classes /, //, /// and IV. All collinear events are contained in class /. 
Other events are distributed according to the number of photon candidates and the apla- 
narity ^apl an- 
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is the angle between photon candidates i and j. Here and in the following sections the photon 
candidates are ordered by energy. The class definitions are shown in Table |[ 

Cuts are applied on the energy sum as well as on the missing transverse and longitudinal 
momenta. The energy sum Es is defined as the sum of the photon candidate energies Ei 
plus the missing longitudinal momentum which could originate from a photon escaping along 
the beam direction. The detailed definitions of these quantities depend on the event class. 
However the cuts summarised in Table ^ namely that the energy sum must be at least 60% 
of the centre-of-mass energy and the missing transverse momentum must be less than 10% 
of the centre-of-mass energy are equivalent. Since the event angle is defined using the two 
highest-energy photons, events are rejected if one of these escapes detection, that is if the 
missing longitudinal momentum is larger than the energy of the second highest-energy photon 
candidate. This cut on the missing longitudinal momentum is not designed to reject background 
but ensures a good signal definition since it prevents ambiguous events from entering the signal 
sample. Distributions of the relevant quantities are shown in Figure |^. 

Class / contains coUinear events which make up about 90% of the total sample. For these 
events it is assumed that missing transverse and longitudinal momenta are negligible, so no 
cuts on these quantities are applied. The energy sum is simply taken to he Eg = Ei + E2. A 
class // event has exactly two observed photon candidates; since the event is acollinear a third 
photon is assumed along the beam direction. Quantities B and -Eiost, equivalent to the missing 
transverse and longitudinal momenta, are calculated assuming three-body kinematics and are 
given by: 

B = ■ (sin^i + sin ^2) I cos - 02)/2]| , (9) 

^lost = v^/ [1 + (sin 01 + sin02)/| sin (^i + ^2)!] ■ (10) 

Consequently Eg = Ei + E2 + -Eiost- For classes /// and IV the missing transverse and longi- 
tudinal momenta, pt and p\, are calculated in the usual way from cluster energies and angles 
and^f = Ei^*+Pi- 



event class 


/ 


II 


III, w 


cut 


energy sum 




E§ 


Ef 


> 0.6v^ 


transverse momentum 




B 


Pt 


< O.lv^ 


longitudinal momentum 




-^lost 


P\ 


<E,,E2 



Table 3: Cuts for the kinematic event selection. The cut variables depend on the class, see text. 
For class / events no cuts on the missing longitudinal and transverse momenta are applied. 



5.4 Neutral event selection 

Events without charged particles produced at the primary interaction point are referred to here 
as neutral events. They are selected using information about the hits in the two main tracking 
chambers CJ and CV. A photon candidate is said to have hits in one of these chambers if more 
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than about 50% of the innermost wires in the sector pointing to the photon have hits. For CV 
the 12 axial layers are taken into account and for CJ the inner 16 layers. 

Events are rejected in either of the following two cases: 

1) double veto: at least two photon candidates have hits in CV. 

2) single veto: any photon candidate has hits in both CV and CJ, unless it is identified as 
a conversion. To be identified as a conversion, information from the SI detector is necessary 
and therefore this condition is restricted to | cos^| < 0.89. A conversion must have exactly two 
tracks reconstructed in a three dimensional cone of half angle 20° around the cluster, with at 
least one of the tracks assigned to the cluster. Additionally, the two SI ladders in the direction 
of the cluster must have no more than two hits. Dead ladders are counted as hits to ensure a 
good background rejection. The two sides of each ladder are counted separately, leading to a 
maximum number of four hits. Events are rejected if any photon with | cos6'| > 0.89 has hits 
in CV and CJ. 



5.5 Summary 

The numbers of events selected and expected after each step of the selection are given in 
Table The preselection is very loose and there are more events observed after the preselection 
than expected from the listed Monte Carlo samples. This excess is mainly due to cosmic- 
ray events and four-fermion events with two observed electrons. The difference between the 
number of observed and expected events becomes much smaller after the rejection of cosmic- 
ray background. After the kinematic cuts the sample consists mainly of Bhabha events and 
the signal efficiency remains at almost 100%. The cut on the longitudinal momentum rejects 
events with a high-energy particle along the beam direction. For Bhabhas these are mainly 
events in which one electron escapes and the second electron and a photon are observed in the 
detector (e7(e) topology). Most of the remaining Bhabha events have both electrons in the 



cut 


data 


SMC 


77(7) 


e"'"e (7) 


e7(e) 


z/i?7(7) 


qq(7) 


T^T (7) 


preselection 


192558 


123751 


5826 


107791 


7280 


105 


398 


2352 


cosmic bkg. 


133099 


122898 


5823 


107697 


7194 


104 


244 


1835 


kinematic cuts 


120515 


119674 


5809 


107048 


6310 


6.6 


130 


370 


longitudinal mom. 


108832 


110082 


5520 


103833 


539 


0.55 


68 


122 


double veto 


6367 


6152 


5505 


68 


515 


0.55 


52 


12 


single veto 


5235 


5261 


5258 


0.38 


1.55 


0.55 


1.10 


0.05 






±12 


±12 


±0.19 


±0.43 


±0.09 


±0.24 


±0.03 



Table 4: The number of events observed in data after the cuts, the signal expectation and the 
most important background sources indicated by their final state are given. The row labelled 
kinematic cuts contains only the cuts on the energy sum and the missing transverse momentum; 
the cut on the missing longitudinal momentum is listed separately. The neutral event selection 
is split up into the double veto and the single veto. For the numbers of events in the final 
selection after the single veto, the statistical error is also given. All Monte Carlo predictions 
are normalised to the integrated luminosity of the data. 
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detector and are easily rejected by the double veto. After the single veto, all background levels 
are almost negligible. Besides the requirement on the longitudinal momentum, most of the 
efficiency loss comes from the single veto, since the double veto rejects only events with two 
converted photons. 



6 Experimental systematic errors 



The efficiency of the selection is studied with a signal Monte Carlo simulation including full 
detector simulation. It is found to be independent of the centre-of-mass energy and therefore 
the average efficiency is used to correct the data. The efficiency is about 97.8% in the barrel 
region and drops to 69% at the edge of the selection at cos6'* = 0.93 as can be seen in Table ||. 
The total efficiency within cos 6'* < 0.93 is 92.6%. Uncertainties in the simulation give rise to 
the systematic errors as discussed below. The errors are given in % relative to the efficiency. 



6.1 Cut values 



To assess the stability of the cuts, the cut values are varied and the difference between data and 
Monte Carlo expectation is assigned as a systematic error. The kinematic cuts on energy sum, 
missing transverse and longitudinal momenta are changed by ±10%. The cuts on the numbers 
of hits allowed in the neutral event selection are varied by one or two hits. The resulting 
systematic errors are at most 0.1% and in total contribute 0.17%. 



6.2 Conversion probability 

A crucial point for this analysis is the correct modelling of the material in the detector since 
this material leads to photon conversions. If a photon converts before the first active detector 
layer, a reliable distinction between a primary produced electron and a photon conversion is not 
possible. Events with such photons must be rejected from the analysis. Events with photons 
which convert later are selected and do not influence the efficiency. 

The probabilities given in this paragraph are conversion probabilities for single photons and 
the given errors are absolute numbers. For | cos^| < 0.89 the Monte Carlo prediction for the 
probability for a single photon to convert before SI is 0.7%. A systematic error of 0.2% is 
assigned to this value since the only material in front of the SI detector is the beam pipe which 
is uniform and therefore well modelled. There is, however, more material in the endcap e.g. from 
SI and CV readout electronics. In the region 0.89 < | cos^| < 0.93 the Monte Carlo prediction 
for the conversion probability in front of CV is 6 - 9%. Studies, for example, of the width 
of the ECAL energy distribution in Bhabha events show no signiflcant disagreement in the 
Monte Carlo description of material in this region, and a 1% error is assigned to the conversion 
probability. For | cos^| > 0.93 signiflcant effects due to unmodelled material are clearly visible 
in these studies and this is the main reason for restricting the two highest-energy photons to 
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jcos^l < 0.93. A 3% systematic error is assigned for additional photons with |cos0| > 0.93. 
The Monte Carlo prediction for the conversion probability in this region is 10 - 15%. 

Most events have two back-to-back photons. These photons hit equivalent kinds of material 
in the left and right side of the detector and consequently their conversion probabilities are 
correlated. For a two-photon event this leads to a systematic error on the efficiency of 0.4% in 
the region | cos6'| < 0.89. For a typical mixture of events at all angles the average systematic 
error on the efficiency is between 0.45% and 2.01%, where the 0.45% is taken to be correlated 
between cos 6** bins. 



6.3 Angular reconstruction 

Another important source of uncertainty is the reconstruction of the photon angle. A systematic 
mis-reconstruction of the angle would lead to distortions in the measured angular distribution. 
Monte Carlo studies show that the angular resolution is approximately 0.2° — 0.3°, which makes 
a full unfolding unnecessary. However, some systematic shifts of up to 0.15° are predicted 
leading to a bias of the measured angular distribution. The Monte Carlo prediction of this bias 
is corrected via the efficiency bin-by-bin and is included in the efficiency given in Table ^ Since 
the angular reconstruction for photons cannot be cross-checked with the data the full size of 
the predicted bias is taken as a systematic error on the efficiency. For most of the central region 
this error is negligible, at other angles the error is between 2% and 3%. For the calculation 
of the total cross-section, these effects are only relevant at the edge of the phase-space of the 
selection, leading to an error of 0.46% corresponding to a shift of 0.15° at cosO* = 0.93. 

While the event angle is mostly affected by shifts that are symmetric on both sides of the 
detector, the acollinearity is changed by shifts that are different on the left and right side. The 
error arising from such shifts is 0.1°. 



6.4 Background 

The expected background is very small. The Monte Carlo expectation for Standard Model 
processes, summarised in Table ^, is less than 4 events. Selected events in which one photon 
has hits in CV are scanned and the Standard Model background is estimated to be less than 
10 events. Background from cosmic-ray events is estimated to be also less than 10 events by 
scanning selected events with at least one track segment in the muon chambers. For both these 
sources half of this estimate is corrected while the other half is taken as a systematic error. 
Since the angular distribution of the background is poorly known, to simplify the correction 
procedure, the shape of the signal distribution is assumed. Thus the background is corrected 
by increasing the efficiency by 0.2% corresponding to 10 events. A systematic error of 0.1% for 
each of the two sources is assigned. 
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6.5 Other errors 

The efficiency loss due to overlaid cosmic-ray events^ and noise hits in the detector is studied 
with randomly triggered beam crossings. This leads to expected signal rejection rates of 0.22% 
due to cosmic-ray events and 0.30% due to noise, corresponding to 12 and 16 events respectively. 
These numbers are cross-checked by scanning events and good agreement is found. The resulting 
efficiency loss of 0.52% is corrected for and a systematic error of 0.05% is assigned for each 
source. 

The error due to Monte Carlo statistics enters only in the determination of the efficiency which 
leads to a very small binomial error of 0.1%. The radiative corrections are determined from 
3 ■ 10^ events, making the statistical error much smaller than the theoretical uncertainty and it 
is therefore neglected. 

The luminosity errors are given in Table |l]. The errors change slightly from energy to energy 
and are dominated by the systematic error. A total common error of 0.23% is assigned in order 
to simplify the error treatment in the log likelihood fits. 

The trigger efficiency of the ECAL is tested with Bhabha events and is found to be 100%. The 
corresponding error is negligible. 

6.6 Summary 

In general the experimental systematic errors for this analysis are small. The largest errors 
arise from the conversion probability, the angular reconstruction and the luminosity. The 
overall systematic error for the total cross-section is 0.77%. For the differential cross-section 
the systematic error as a function of cos 9* is given in Table |^. For most of the central region 
the total systematic error is 0.56%, which is composed of 0.45% from the conversion probability, 
0.23% from the luminosity, 0.17% from the cut values and small errors from the background 
and the efficiency loss. All systematic errors are independent of the centre-of-mass energy and 
fully correlated between energies. 



7 Results 

7.1 Cross-sections 

The number of observed events in each class is given in Table |^ and compared to the number of 
events expected from Monte Carlo simulation. There is a good agreement for collinear events 
(class /), but there are fewer acoUinear events observed than expected. To judge whether this 
is an indication for physics beyond the Standard Model one has to take into consideration that 
the main Monte Carlo generator includes only 0{a^) terms. There is no Monte Carlo program 
available to calculate events with four photons in which one travels along the beam direction, 

•^A cosmic-ray particle in the detector which coincides with a signal event. 
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class 


/ 


II 


III 




IV 




> 2 


2 


3 


3 


> 4 


observed 


4747 


394 


71 


20 


3 


expected 


4730 


435 


94 




5.3 



Table 5: The number of observed and expected events in the different classes. The Monte Carlo 
prediction for classes / - /// is from RADCOR. The expectation for class IV with four photons 
is calculated using FGAM. There is no Monte Carlo program available to generate events with 
three observed photons plus one escaping along the beam direction, hence no prediction for 
class IV events with three observed photons is given. 



so no expectation is given for class IV events with three observed photons. Events with four 
observed photons are calculated using FGAM [p^ . 

The most important discriminant for the classes is the acoUinearity, the distribution of which 
is shown in Figure |^. The typical resolution is 0.35°. Discrepancies between data and Monte 
Carlo simulation at ^acoi < 1° may occur for two reasons. One is that Monte Carlo events with 
a soft third photon below some energy cut-off are generated with two exactly back-to-back 
photons. This affects the acoUinearity distribution at ^acoi < 1° as can be seen in Figure |^. 
The other reason is the possible systematic shift of the reconstructed angle of about 0.1°. For 
larger values of ^acoi differences are most likely due to missing higher order effects. The largest 
significance of a discrepancy between data and Monte Carlo simulation occurs for a cut of 
^acoi > 3°. In this case 1141 events are predicted with 1004 observed. Including a systematic 
error of 16 events due to the mismodelling of the angle, an excess in the Monte Carlo prediction 
of 3.6 standard deviations is observed. This discrepancy may be due to higher order effects, 
which can be large. 

The measured differential cross-sections are plotted in Figure § for six energy ranges. Infor- 
mation for all eight energies can be found in Table |^. Integrating these distributions leads to 
the total cross-sections which are shown as a function of ^/s in Figure ^ and are summarised 
in Table 0. The cross-sections are in very good agreement with the QED expectation. On 
average cXobs/cQED = 0.999 ± 0.014 ± 0.008, where the first error is statistical and the second 
is systematic. This average has a x^/dof = 12.7/7 and includes the correlation of systematic 
errors. There is an additional error from theory of 0.01. 

7.2 Tests of alternative models 

Binned log likelihood fits for the alternative models are performed on the measured differential 
cross-sections taking systematic errors and their correlations into account. Where possible, the 
fit parameters, given in Table ||, are chosen such that the resulting probability distribution 
is approximately Gaussian. One-sided 95% confidence level limits are obtained for the fit 
parameters by renormalising the likelihood function to the physically allowed region in the 
same way as in 0. The quoted limits on the model parameters correspond to these limits on 
the fit parameters. For excited electrons two types of fits are performed. Limits are determined 
on the coupling constant /^/A of an excited electron to an electron and a photon for various 
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183 


189 


192 


196 


200 


202 


205 


207 


efficiency 


r ad. cor. 

/V 


cosy -Dm 


events 


0.00-0.05 


8 


33 


3 


10 


15 


8 


17 


23 


0.9515 


0.0250 


1.0689 


0.05-0.10 


21 


37 


5 


18 


10 


6 


14 


20 


0.9769 


0.0056 


1.0666 


0.10-0.15 


15 


34 


4 


14 


15 


4 


13 


15 


0.9769 


0.0056 


1.0644 


0.15-0.20 


10 


40 


7 


15 


12 


3 


11 


22 


0.9769 


0.0056 


1.0621 


0.20-0.25 


13 


40 


2 


9 


13 


5 


14 


23 


0.9769 


0.0056 


1.0599 


0.25-0.30 


14 


43 


10 


11 


13 


2 


8 


27 


0.9769 


0.0056 


1.0577 


0.30-0.35 


20 


47 


10 


14 


20 


7 


29 


26 


0.9755 


0.0056 


1.0556 


0.35-0.40 


22 


41 


15 


22 


21 


17 


12 


23 


0.9753 


0.0056 


1.0534 


0.40-0.45 


21 


38 


12 


19 


11 


9 


18 


18 


0.9753 


0.0056 


1.0512 


0.45-0.50 


13 


49 


8 


22 


15 


6 


17 


31 


0.9753 


0.0056 


1.0492 


0.50-0.55 


21 


61 


7 


23 


24 


10 


22 


23 


0.9742 


0.0056 


1.0471 


0.55-0.60 


23 


72 


6 


25 


25 


18 


21 


37 


0.9742 


0.0056 


1.0451 


0.60-0.65 


25 


64 


14 


38 


24 


19 


21 


48 


0.9742 


0.0056 


1.0433 


0.65-0.70 


24 


96 


18 


35 


33 


17 


37 


68 


0.9762 


0.0056 


1.0415 


0.70—0.75 


36 


104 


18 


42 


34 


23 


50 


75 


0.9906 


0.0208 


1 A OA A 

1.0399 


0.75-0.80 


48 


131 


23 


59 


53 


19 


47 


81 


0.9824 


0.0208 


1.0387 


0.80-0.85 


66 


176 


33 


71 


72 


35 


71 


89 


0.9647 


0.0208 


1.0380 


0.85-0.90 


78 


255 


32 


102 


75 


48 


84 


145 


0.8976 


0.0220 


1.0387 


0.90-0.93 


60 


170 


31 


67 


69 


25 


60 


97 


0.6937 


0.0286 


1.0412 



Table 6: Measured angular distributions. For each centre-of-mass energy and cos^* bin the 
number of observed events is given. The last columns give the efficiency, its relative systematic 
error and the radiative correction which has an error of 0.01. The efficiency includes the 
corrections due to the background and the rejection because of noise and overlaid cosmic-ray 
events. For the systematic error on the efficiency the common contribution of 0.56% is taken 
as correlated between bins. 





Born cross-section [pb] 


[GeV] 


observed 


QED 


182.68 


10.05±0.43±0.08 


9.32 


188.63 


8.79±0.23±0.07 


8.74 


191.59 


9.24±0.58±0.07 


8.47 


195.53 


8.43±0.34±0.07 


8.13 


199.52 


7.39±0.31±0.06 


7.81 


201.63 


7.88±0.47±0.06 


7.65 


204.88 


7.40±0.31±0.06 


7.42 


206.56 


6.78±0.23±0.05 


7.29 



Table 7: The measured and predicted total cross-sections at Born level within the angular range 
of cos^* < 0.93. The measured values are shown with their statistical and systematic errors. 
The additional uncertainty on the theoretical prediction is about 1%. 
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masses Mg*. The resulting limits shown in Figure |^ do not depend on a special choice of A. 
A limit on the mass Mg* for a fixed coupling constant = 1 is also determined. In this case 
the scale A is fixed to Mg* and no fit parameter can be found to give a Gaussian probability 
distribution. Therefore the likelihood distribution is shown in Figure 0. To make the Standard 
Model value at Mg. = oo visible the likelihood is plotted as a function of M^^. The limit of 
Mg* > 245 GeV is obtained at ALogL = -lnC + In^^ax = 1-92. 

All fit results are summarised in Table ^. The limits on the model parameters are stronger by 
10% to 20% than our previously published limits 0. Limits on the excited electron coupling 



fj/A are much weaker than the limits obtained from direct searches |^0| but extend beyond 



the kinematic limit for excited electron production as can be seen in Figure |[ As a comparison 
limits from excited electron production via electron-7 fusion determined by HI are also 
shown. Similar results are available from ZEUS |H2 . 



Fit parameter 


Fit result 


95% CL Limit [GeV] 




(-40.3+^^:^) TeV-^ 


A+ > 371 
A_ > 314 




[-3.56tijj) TeV-6 


A' > 800 




(0.926^S:S) TeV- 


A = +1: Ms > 805 
A = -1: Ms > 956 


Me.(^ = l; A = Mg*) 


see Figure 


Mg* > 245 


(/^/A)4(Mg. = 200 GeV) 


(-215^-) TeV- 


f^/A < 4.11 TeV"^ 



Table 8: Fit results and limits at 95% confidence level obtained from binned log likelihood fits 
to the differential cross-sections. The model parameters are defined in Section |. 



7.3 Resonance production 

Of the 71 selected class III events, 64 have all three photons within |cos^| < 0.93. These 
events are used to search for a photonically decaying resonance X which is produced in asso- 
ciation with a photon (e^e~ — > 7X , X 77). Each event is a candidate for three different 
masses corresponding to the pairing of the photons. For planar three photon events three-body 
kinematics can be used to calculate the energies Ek of the photons: 

Ek oc sin Oij ; Ei + E2 + E^ = yfs , (11) 

where Ek is the energy of one photon and Uij is the angle between the other two photons. This 
leads to a resolution on the invariant mass of photon pairs of about 0.5 GeV. Calculating the 
mass from the observed cluster energies would result instead in a resolution of about 3%. 

For the signal it is assumed that the resonance X has isotropic production and decay distribu- 
tions. The QED background on the other hand is peaked in the forward direction. Since the 
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[GeV] 


Mass [GeV] 


[GeV] 


Mass [GeV] 


181.73 


51.7 


84.1 


152.6 


195.47 


100.3 


108.4 


128.1 


182.41 


47.7 


83.5 


155.0 


199.55 


56.2 


126.9 


143.3 


182.72 


55.7 


73.8 


157.6 


199.56 


43.1 


114.3 


157.8 


182.75 


58.7 


49.5 


165.8 


199.57 


81.3 


48.3 


175.7 


182.75 


42.7 


68.4 


164.0 


195.59 


56.0 


106.5 


154.2 


182.86 


60.1 


118.2 


125.9 


195.59 


80.3 


82.9 


158.0 


182.89 


89.0 


104.2 


121.1 


199.41 


90.0 


96.8 


149.3 


182.70 


51.3 


103.2 


141.8 


199.45 


78.3 


121.1 


137.8 


182.70 


67.6 


116.3 


123.6 


199.52 


99.3 


116.3 


128.1 


182.70 


69.5 


84.1 


146.5 


201.58 


40.8 


127.3 


150.9 


188.58 


60.4 


50.7 


171.3 


201.63 


81.8 


118.9 


140.8 


188.59 


43.5 


88.8 


160.6 


202.55 


49.3 


99.4 


169.5 


188.59 


79.9 


69.0 


156.2 


203.66 


83.1 


78.0 


168.8 


188.59 


102.2 


78.8 


137.5 


204.65 


26.2 


127.3 


158.1 


188.60 


33.6 


94.8 


159.5 


204.73 


51.4 


86.8 


178.2 


188.61 


76.3 


102.7 


138.5 


205.04 


54.9 


64.9 


186.6 


188.62 


49.5 


73.7 


166.4 


205.11 


60.5 


78.2 


179.7 


188.63 


66.5 


45.9 


170.4 


205.13 


65.0 


84.6 


175.2 


188.63 


51.1 


76.0 


164.9 


205.14 


62.5 


93.9 


171.4 


188.64 


105.7 


97.8 


121.8 


205.35 


37.7 


137.0 


148.3 


188.65 


37.9 


98.7 


156.3 


205.81 


73.1 


120.4 


150.0 


188.66 


115.3 


99.0 


111.8 


206.03 


62.4 


84.4 


177.3 


188.66 


35.7 


94.4 


159.4 


206.17 


63.1 


64.2 


185.5 


188.66 


67.6 


65.8 


163.4 


206.17 


58.4 


103.1 


168.7 


188.67 


93.0 


77.9 


144.5 


206.39 


54.2 


92.4 


176.4 


188.67 


35.1 


96.4 


158.3 


206.55 


66.0 


116.7 


157.2 


188.67 


54.7 


86.3 


158.6 


206.56 


93.8 


104.1 


151.8 


188.67 


65.4 


92.8 


150.7 


206.57 


83.3 


65.0 


177.5 


189.01 


34.4 


99.3 


157.1 


206.58 


105.5 


125.5 


125.7 


191.52 


41.2 


84.4 


166.9 


206.61 


69.6 


108.7 


161.4 


191.63 


43.2 


116.9 


145.6 


207.96 


69.9 


84.9 


176.5 


195.47 


61.5 


83.3 


165.8 


208.03 


43.1 


112.1 


169.9 



Table 9: The three invariant masses for planar three photon events. The events are grouped 
according to the energy ranges given in Table |l]. 
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QED Monte Carlo generator appears to overestimate the production of three-photon events the 
background is scaled to the number of observed events. 



The efficiency for resonance production can be separated into two parts. The first part accounts 



for the restricted phase-space of the selection summarised in Equation 12 



I cos I < 0.93, z = 1,2,3 ; ^acoi > 10° . (12) 

This efficiency depends strongly on the mass, Mx, of the resonance. At a centre-of-mass energy 
of 207 GeV it is between 64% and 72% for masses in the range 70 GeV to 170 GeV. The mass 
range is limited by the acoUinearity cut. The second part is the reconstruction efficiency for 
X7 events within the phase-space defined by Equation this depends only weakly on Mx 
and the centre-of-mass energy. It is 96.3% which is larger than the efficiency for QED events 
because of the different angular distribution. 

The observed distribution of the invariant mass of photon pairs is shown in Figure |^ together 
with the distribution expected for a signal at Mx = 131.5 GeV. A list of the three invariant 
masses per event can be found in Table ^. Limits at 95% confidence level on the product 
of X7 production cross-section and branching ratio are derived using a method of fractional 
event counting |^|, where the weight of an event is determined according to the resolution 



and the difference between hypothetical and reconstructed masses. The intrinsic width of the 
resonance X is assumed to be negligible. The limits are calculated assuming that the cross- 
section is independent of i/i. If the production cross-section depends on ^/s, the hmit can be 
used with respect to the luminosity weighted centre-of mass energy of (y/s) =196.6 GeV. For 
most production cross-sections (e.g. a oc 1/s or oc In s) the difference from the correctly scaled 
limit is less than 0.8% for masses in the range 40 GeV to 170 GeV. For a threshold behaviour 
like cr oc (1 — [|TH[ this holds only up to Mx < 145 GeV. The limits are about a factor two 
better than our previously published results 0. 



8 Conclusion 



Total and differential cross-sections for the process e^e~ 77(7) have been measured at high 
energies with a large data set. The data are in good agreement with the QED expectation, 
though some discrepancies in the acoUinearity distribution are observed, which are attributed 
to missing higher order effects in the 0{a^) model calculations. Strong constraints on models 
predicting deviations from QED are obtained which are summarised in Table |[ Lower limits 
at 95% confidence level are placed on cut-off parameters A± of about 340 GeV and on the 
scale of gravity in extra dimensions of about 880 GeV. Excited electrons must be heavier than 
245 GeV if the relative strength of the e*e7 vertex is /7/A = 1/Mc*. Limits on this coupling 
are placed for a range of excited electron masses. In the mass spectrum of photon pairs no 
indication for a narrow resonance X is found leading to limits on cross-section times branching 
ratio for X7 production with X ^ 77 of about 0.02 pb assuming isotropic production and 
decay distributions. 
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Run: 16249, Event: 18045 

Date: 2 November 2000, Time: 07:57:38 
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Side view - plane of Th 




Phi=383. 





ll;^ I 1 \ \ S / de V i ew - plane pe rp . 

X - y V i ew 

Figure 1: Last OPAL e"'"e~ — >■ 77(7) event, taken three minutes before the final shut-down of 
LEP. Two high-energy clusters are detected in the electromagnetic calorimeter. This is a class 
// event with acollinearity ^acoi = 17°. One photon has converted between CV and CJ; the two 
corresponding tracks are visible in the tracking chambers. 
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OPAL 




Figure 2: Distributions used in the kinematic event selection. The energy sum for events of 
all classes is shown for barrel and endcap regions separately in plots a) cos^* < 0.81 and b) 
cos 6* > 0.81. For events that are not in class / plot c) shows the transverse momentum and d) 
the longitudinal momentum. The points are the OPAL data after application of all cuts except 
that on the quantity which is plotted. The histograms show the Monte Carlo expectation. The 
arrows in a), b) and c) show the positions of the cuts. For the longitudinal momentum there 
is no cut at a specific value, events are selected if pi < £'1,2 as indicated by the Monte Carlo 
distribution shown as the shaded histogram. Background comes mainly from cosmic-ray events. 
In plot d) some Bhabha background is visible. 
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^acol [deg] 

Figure 3: AcoUinearity ^acoi for all selected events. Because of the non-uniform bin size the 
entries are normalised to events/l°. The measured distribution is compared to the 0{a^) Monte 
Carlo prediction including full detector simulation. The Monte Carlo distribution at generator 
level is also shown. The resolution is about 0.35° and the possible systematic bias on the angle 
is around 0.1°. 
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Figure 4: The measured differential cross-section at the Born level for the process e^e~ 
77(7) for six ranges of centre-of-mass energy. The points show the number of observed events 
corrected for efficiency and radiative effects. The solid curve corresponds to the Born-level 
QED prediction. The dashed lines represent the 95% confidence level interval from the fit to 
the function given in Equation |^. 
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Figure 5: Total cross-section at the Born level for the process e"'"e~ 77(7) with cos6* < 0.93. 
The curve corresponds to the Born-level QED expectation. The data are corrected for efficiency 
loss and higher-order effects. The errors shown for the measurements are statistical only. The 
systematic errors are correlated between energies and are plotted as a band around the QED 
expectation. The inner band represents the experimental error and the outer band includes the 
theoretical error. 
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Figure 6: Upper limit at 95% confidence level on the coupling /-y/A as a function of the mass 
Me* of an excited electron. 
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Figure 7: Likelihood difference as a function of M~t^ for fixed fy = l and A = Mc* . The limit 
obtained at ALogL = 1.92 is M^* > 245 GeV. 
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Figure 8: Comparison of limits on the coupling of excited electrons. OPAL results are given 
from this analysis and the two channels of the direct search [4C] together with the results from 
HI pi[|. The shaded regions are excluded at 95% confidence level in each case. 
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Figure 9: Results of a search for resonance production in class /// events, a) shows the invari- 
ant mass of photon pairs for data (points) and for the e"'"e~ 777 Monte Carlo expectation 
(histogram) scaled to the number of observed events. There are three photon pairs per event. 
A hypothetical signal with (7x7 x BR = 0.02 pb and a mass Mx = 131.5 GeV is superim- 
posed (hatched histogram). The binning is chosen to match the expected mass resolution, b) 
shows the upper limit at 95% confidence level for the product of production cross-section and 
branching ratio for the process e+e" X7, X ^ 77 as a function of the mass of the resonance 
X. 
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